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ABSTRACT 

When VLBI observations are used to determine the position or motion of 
a radio source relative to reference sources nearby on the sky, the astrometric 
information is usually obtained via: (i) phase-referenced maps; or (ii) paramet- 
ric model fits to measured fringe phases or multiband delays. In this paper we 
describe a "merged" analysis technique which combines some of the most impor- 
tant advantages of these other two approaches. In particular, our merged tech- 
nique combines the superior model-correction capabilities of parametric model 
fits with the ability of phase-referenced maps to yield astrometric measurements 
of sources that are too weak to be used in parametric model fits. We compare the 
results from this merged technique with the results from phase-referenced maps 
and from parametric model fits in the analysis of astrometric VLBI observations 
of the radio-bright star IM Pegasi (HR 8703) and the radio source B2252+172 
nearby on the sky. In these studies we use central-core components of radio 
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sources 3C 454.3 and B2250+194 as our positional references. We obtain astro- 
metric results for IM Peg with our merged technique even when the source is too 
weak to be used in parametric model fits, and we find that our merged technique 
yields superior astrometric results to the phase-referenced mapping technique. 
We used our merged technique to estimate the proper motion and other astro- 
metric parameters of IM Peg in support of the NASA/Stanford Gravity Probe B 
mission. 

Subject headings: astrometry — binaries: close — radio continuum: stars - 
stars: activity — stars: imaging — stars: individual (IM Pegasi) — techniques: 
interferometric 



1. Introduction 



Very-long-baseline interferometry (VLBI) provides the most accurate astrometric mea- 
surements of celestial objects currently attainable. VLBI astrometry has been used, among 
other applications, to: define the most nearly inertial referen ce frame available for the posi- 
tions of celestial objects (e.g-. lMa et al.lll998t iFey et al.l 120091 ): tie an iner tial (extragalactic) 



reference frame to a p l anetary ephemeris y ia observations of pulsars (e.g., iBartel et al.lll985 



Rodin fe Sekido 



Earth (e.g., iRyan et al 



2002 



1993 



Dodson et al.ll2003l); characterize m otions and other properties of the 



Mathews &: Shapiro! Il995l ); study posit ions and motions of 



maser spots in galaxies as a means to estimate black-hole masses fe.g.. Ilshihara et al.ll2001 



Kondratko et al.l 120041 ) and distances to other galaxies (jHerrnstein et al. 



general relativity via mea surements of solar gravitational deflection (e.g.. ICounselman et al. 



1999): and test 



1974 : IShapiro et al.ll2004j ). In this paper we describe the analysis of astrometric VLBI obser- 
vations of the RS CVn binary star IM Pegasi (IM Peg; HR 8703), the radio-bright star which 
served as the "guide star," and hence as the positional reference, for the NASA/Stanford 
Gravity Probe B (GP-B) experiment. 

GP-B was designed to m easure the geodetic and frame- dragging effects predicted by 
general relativity (see Paper I. IShapiro et al.ll201ll ). as manifested in secular changes in the 
spin-axis orientations of four gyroscopes placed within a spacecraft in a low-altitude, polar 
orbit about the Earth. The spacecraft also had an on-board telescope equipped with a 
tracking system designed to keep the guide star, IM Peg, at the center of the telescope field 
of view, and thereby provide a directional reference: the relativistic effects were estimated 
by measuring the drift rates in the spin-axis directions of the four gyroscopes relative to the 
direction to IM Peg. Thus, if any motions of IM Peg on the sky relative to an inertial frame 
were not accounted for sufficiently accurately, they could map directly into the gyro drift- 



-3 - 



rate signals and thereby corrupt the relativistic measurements. (Ideally the GP-B guide star 
would have been a distant quasar, which could be treated as effectively motionless on the sky, 
rather than a star in our galaxy, but the on-board telescope could only track objects brighter 
than about 6th magnitude, so a relatively bright star was the only possible choice. A list of 
criteria that had to be met by th e guide s tar, and a description of the selection process that 



led to the choice of IM Peg, is in lPaper J .) Our VLBI observations of IM Peg were therefore 
undertaken to determine the motions of the star relative to an inertial frame so that the 
relativistic effects included in the gyro measurements could be properly separated. The goal 
of our VLBI observations was to determine, relative to an inertial frame, the proper motion 
of IM Peg with standard error 0.14 mas yr _1 or less in each of the north-south and east- 
west directions. This accuracy goal was based upon an error projection for the full GP-B 
experiment prior to launch in April 2004 and the desire at the time for the uncertainty in the 
proper motion of IM Peg not to constitute a significant source of error for the experiment. 

In this paper we present a new approach for astrometric VLBI data analysis that we 
developed to meet the accuracy goal for GP-B. We then compare the results from this ap- 
proach with the results from two well-established VLBI analysis techniques: phase-referenced 
mapping (PRM) and "phase connection" followed by parametric model fitting (PMF). Our 
new approach is a combination of these two standard approaches that provides many of the 
benefits of each. All three approaches use differential VLBI astrometry, in which the target 
source (i.e., the source of astrometric interest) is observed alternately with at least one ex- 
tragalactic, compact reference source nearby on the sky to provide model corrections in the 
analysis of the target-source data. The way these corrections are applied differs among the 
approaches, but in all cases the astrometric measurements of the target source are relative 
to the position(s) of the reference source(s). 

In § H] we briefly describe the VLBI observables used in our analysis. We then discuss 
some basic aspects of the PRM and PMF analysis techniques in § Eland §SJ respectively. In 
§[5]we describe our new "merged" analysis technique and then our implementation of it. We 
compare the results from the three techniques in § [61 and offer corresponding conclusions in 
§0 



2. VLBI Observables: Fringe Phase, Fringe Rate, and Multiband Delay 

In VLBI observations of a compact source, the total fringe phase (or visibility phase, 
or interferometric phase, or simply "phase"), (p(u,t), is the phase at frequency u and time t 
that is associated with the difference in arrival times of signals received from the source at 
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two antennas in a VLBI array. We describe 4>(oj, t) as: 

(f)(u, t)=U) [r gcom (t) + T inst (w, £) + T atm (w, t) + T stTUC (u, t) + T noisc (u, t)] + 27rA^(w, t) (1) 

where r geom (t), the "geometric delay," is the difference in the signal arrival times in vacuum 
at the two antennas; r inst (a;,t) represents the difference in the instrumental delays (includ- 
ing clock behavior) at the two antenna sites; r atm (u;,£) represents the difference in signal 
propagation times to the two antennas due to all atmospheric effects, including the contri- 
butions of the hydrostatic (or "dry") atmospheric constituents, atmospheric water vapor, 
and the ionosphere; r stTUC (uj, t) is the delay contribution from source structure, i.e., from the 
non-pointlike brightness distribution of the source; T QO i se (u,t) represents the (thermal) noise 
contribution to the phase measurement; and N(u, t) represents the integer number of "ambi- 
guities," or "phase wraps," included in the measurement. Information about the number of 
phase wraps associated with a fringe phase is not inherent to a measurement of fringe phase, 
and in general a measured fringe phase can be defined such that — % < meas (u;,t) < it, or 
< <fimeas(u,t) < 2%. There is a fringe phase associated with each antenna pair, or "base- 
line," within the VLBI array. The "fringe rate" is the partial derivative of the fringe phase 
with respect to time. The "multiband delay" is a measured approximation of the group 
delay, which is the partial derivative of the fringe phase with respect to frequency. 

To lowest order (i.e., neglecting, e.g., relativistic effects), the geometric delay is given 

by: 

Tgeom(t) = - [B(t) • §(t)] (2) 

where c is the speed of light in vacuum, B(£) is the 3-dimensional vector between the two 
antennas of the baseline, and s(t) is the unit vector in the direction of the observed source. 
All of the astrometric information in the measured fringe phase is contained in T geom (t). 
(While we do not show the relativistic contributions to r gcom (t) in Equation [2], we do include 
those contributions throughout our analyses. We also account for Earth motions relative to 
the solar-system barycenter during the time between signal arrival at the two antennas.) 

The model we used for T atm (u;,t) can be more explicitly described by: 

Tatm A (w, t, e, s) = r zcn _ dryA (t)m dryA (t, t) + T zcn _ wetA (t)m wetA (t, e) + r ionA (u, t, s) (3) 

where the 11 A" subscript refers to site A; r zcn _ dryA (t) is the propagation delay through the 
atmosphere at zenith under the assumption that the atmosphere is in hydrostatic equilibrium; 
T Z en_wct A (£) is the addition al propagation de lay through the atmosphere at zenith due to 



tropospheric water vapor (jDavis et al.lll985l ): md ryA (t,e) and m wctA (£, e) are, respectively, 



the "mapping functions" that project (i.e., scale) the "dry" and "wet" delays at zenith to 
the line-of-site elevation, e, of the observed source; and r ionA (a;, t, s) denotes the line-of-site 
contribution of the ionosphere to phase delays in the direction of the observed source. 
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3. Phase-Referenced Mapping (PRM) 



The basic idea behind the PRM technique is to use the data from the observed reference 
sources to account for otherwise unmodeled fluctuations in instrumental or atmospheric de- 
lays, as well as for other model errors. Specifically, the "residual" fringe phases and rates, 
i.e., the differences between the measured and a priori model values of these quantities, are 
obtained for the reference sources and then temporally and sometimes also spatially interpo- 
lated to the observation time and sky position of the t arget source to estimate the effects of 



model errors on the target- s ource observables (see, e.g . . IShapiro et al. 



19831 ; iLestrade et al.l Il990l ; iBeasley fc Conway! Il995l ; iFomalontl 120051 ) . Typically the cycle 



1979: Gorenstein et al. 



time over which the reference sources and target source are observed is relatively short, from 
several seconds to several minutes. To the extent that reference-source structure is properly 
accounted for and the reference-source residuals are properly interpolated to the observation 
time and sky position of the target source, the remaining residual components in the target- 
source data can be attributed primarily to measurement noise, unmodeled target-source 
structure, and a position offset of the target source relative to the a priori model position. 
Perhaps the most important feature of the PRM technique is that the target-source data 
can be coherently integrated over the entire span of an observing session. Thus, even very 
weak sources with flux densities well under 1 mJy can be detected and imaged with VLBI 
via this technique. 

Henceforth in this paper, a "scan" refers to a single continuous observation of a particular 
source. Fringe-phase and fringe-rate measurements from successive reference-source scans 
can provide proper model fringe-phase adjustments for intervening target-source data only 
if the change in the model errors between the scans is <C 2n, so that the number of 2-k 
phase wraps between successive reference-source scans can be accurately tracked. For that 
reason, the cycle time over which the reference and target sources are observed should be 
as short as possible. On the other hand, the reference-source scans must be long enough 
to provide reliable fringe-phase measurements from single scans, and the scan times for the 
target source should be long enough to ensure sufficient coherent integration time over the 
course of the observing session to produce an image with acceptably high signal-to-noise 
ratio (SNR). These trade-offs must be balanced to determine the "optimum" scan and cycle 
times in the observation schedule. 

Another beneficial feature of the PRM approach is that it is relatively quick and efficient 
for obtaining high-accuracy astrometry, in large part due to the software packages now readily 
available. In our analyses we used almost exclusively the Astronomical Image Processing 
System (AIPS) provided by the National Radio Astronom y Observatory (N RAO) to produce 
our phase- referenced maps. We followed the guidelines of Diamond! (119951 ). Specifically, we: 
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1. Incorporated the instrumental phase and amplitude calibrations routinely provided 
in log files from individual antennas. We also unweighted faulty data based upon 
information in these files and from operator and correlator reports. 

2. Incorporated an additional constant (over time) phase adjustment to the calibration 
for each antenna based upon "fringe fitting" (via AIPS task FRING) of data from a 
selected reference scan0 Sometimes we used multiple reference scans to provide such 
phase calibration adjustments for all antennas. 

3. Ran a "global fringe fit" (via AIPS task FRING) on all reference-source data while 
applying the phase and amplitude calibrations obtained so far. 

4. Ran AIPS task IMAGR, with the fits from Step 3 as calibration, to generate "self-cal" 



maps of the reference sources (see Paper II. lRansom et al.ll2011al ). We then used these 



self-cal maps to refine the amplitude and phase calibration as well as the data flagging. 
In general we repeated Steps 3 and 4 several times. 

5. Used AIPS task BPASS to further refine the phase and amplitude calibrations as a 
function of frequency. 

6. Used AIPS task IMAGR with the final calibrations from Steps 1 through 5 to obtain 
a phase-referenced map of the target source. 

Unfortunately, the use of the PRM technique for astrometry also poses some challenges, 
especially for campaigns such as ours in which the observations are made over many years. 
For example, one must take great care to assure that model components such as site posi- 
tions, antenna axis offsets, and Earth orientation parameters (EOPs, i.e., X- and Y- pole 



posit i ons, U T1— UTC, and nutation angles in longitude and obliquity; see, e.g., ISeidelmann 



19821 . 119921 ) are modeled consistently and correctly throughout all observation sessions. (For 
example, if updated values of axis offsets for antennas are used in processing data from 
later observing sessions, then the new values will have to be incorporated into a re-analysis 
of data from earlier observing sessions to avoid possible systematic errors.) Furthermore, 



1 The use of suc h constant phase adjustments based on data from selected scans is referred to as "manual" 
phase calibration ( Diamond! 19951 ). Even when measured phase calibrations were available in log files to 
account for instrumental effects, we found that we obtained better calibration across our observed bandwidth 
when we used manual phase calibrations after the application of our measured phase calibrations. This finding 
was based upon an assessment of the phase scatter across the spanned bandwidth of data from some of our 
f 997 and 1998 observations. As a result, we used manual phase calibrations for all data after we applied any 
available measured phase calibrations. 
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post-processing adjustments to model parameters can be quite cumbersome in some analysis 
packages, including AIPS, especially if the adjustments have to be made to a large num- 
ber of experiments. In addition, without the use of multiple reference sources and special 
interpolation routines to handle the spatial (in particular, elevation-angle) dependencies of 
atmospheric delays, a simple temporal interpolation of reference- source fringe phases intro- 
duces model-correction errors that scale roughly with the angular separation between the 
reference source and the target source. Thus the basic PRM technique commonly loses 
viability for an angular separation between reference and target sources larger than a few 
degrees. Finally, the PRM technique offers no inherently good way to assess the effects of 
systematic errors (such as those due to inaccurate modeling of the atmosphere). Reliable 
estimates of t he true accuracy of the astrometric results can therefore be difficult to obtain 
(although see iPradel et al.l 120061 for a comprehensive assessment specifically for the VLBA 
and EVN arrays). Our development and use of the "merged" technique described in §E]of 
this paper was motivated in part by our desire to overcome these drawbacks. 



4. Parametric Model Fitting (PMF) 

The PMF technique for analysis of VLBI data is, in all implementations with which 
we are familiar, essentially the method of "weighted-least-squares" estimation. For each 
reference-source and target-source scan, the total measured fringe phase and fringe rate are 
estimated for each baseline at an epoch near the center of the scan period. Any 2n "phase 
jumps" between successive phase measurements are then resolved (or flagged as unresolvable) 
for each source and baseline. The overall integer number of 2ir differences between the phases 
from different sources is also resolved. Finally, the complete collection of "phase-connected" 
data is fit to a model in which corrections to a wide range of parameters — including in- 
strumental delays, propagation delays through the atmosphere, and source positions — are 
simultaneously estimated. 

To avoid the often onerous task of determining the change in the number of 2n phase 



wrap s between successive scans prior to model fitting, multiband delays (e.g., Park et al. 



19851 ) are commonly used in place of fringe phases as the principal observables with the 
PMF technique. The use of multiband delays comes at a significant cost in statistical (i.e., 
SNR-derived) measurement accuracy relative to the use of fringe phases, but the use of only 
bright radio sources as well as observations over wide spectral bandwidths can reduce this 
cost. The PMF approach with multiband delays is used by the geodetic VLBI community 
in virtually all of their data analyses, including their studies to define the Int ernational 



Celestial Reference Frame and track apparent motions of extragalactic sources (IMa et al. 
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19981 ; iFey et al.ll2009l ). In general, the radio sources chosen for such studies are the brightest 



and most compact known. 

The PMF technique has several advantages over the PRM technique. For example, 
since total measured fringe phases and rates (rather than residual fringe phases and rates) 
are processed in the analysis, one can make improvements or changes to a priori models used 
with multi-year data sets much more easily. Furthermore, the use of data from multiple 
reference sources to account for spatial (e.g., elevation-angle-dependent) as well as temporal 
phase variations is straightforward and readily implemented with existing software packages. 
In addition, because the PMF technique is a method of weighted-least-squares estimation 
of parameters that are treated as Gaussian random variables, one can readily obtain a 
measure of the sensitivity of one estimated parameter (such as a source declination) to the 
variation in another (such as the propagation delay through a site atmosphere) by evaluating 
the covariance matrix obtained with the parameter estimates. The effects on estimated 
parameters of other model changes are also straightforward to evaluate, so various sensitivity 
studies to assess systematic errors can be readily implemented. 

Unfortunately, the PMF technique also has its limitations and problems. Unlike the 
PRM technique, the PMF technique can make use of only those data for which fringe phases 
for a baseline are reliably detected in a single scan. Thus the technique is only well suited 
for relatively bright sources. Furthermore, the PMF technique offers no direct mechanism 
to account for sou rce structure. One can assess fringe amplitudes and "closure" phases 



( IRogers et al.lll974j ) to identify the presence of significant structure, but the model corrections 
for structure have to be generated in a separate process. Finally, when accuracy requirements 
dictate that the PMF technique be used with fringe phases rather than multiband delays, the 
integral number of 2tt phase wraps between neighboring scans must be determined, which 
can be a very labor-intensive process and sometimes is not even possible. 



5. New Merged Analysis (MA) Technique 
5.1. Motivation 

We desired the most accurate astrometry that we could attain from our VLBI observa- 
tions of IM Peg. Unfortunately, neither the PRM nor the PMF analysis technique adequately 
met our needs. 

The PRM technique, as implemented with AIPS, relied upon insufficiently accurate 
models for our astro metric demands, in particular for r atm (o;, t) and also for EOPs and other 
parameters. These models were difficult to reliably correct properly within AIPS. Also, the 
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use of correct, consistent values for some model parameters, such as antenna coordinates 
and axis offsets, would have been somewhat burdensome to implement and assure over the 
8.5 year span of our VLBI observations (which are described further in § 15.3. lj) . In addition, 
we wanted to make full use of observations of a second reference source to model the effects 
of atmospheric gradients on our measured fringe phases, and we had no satisfactory software 
tools with this capacity available to us with the PRM technique^] Finally our use of the PRM 
technique made difficult a robust assessment of several possible sources of systematic error 
in our astrometric results, including, for example, our sensitivity to errors in our a priori 
EOP values or in the various components of our atmospheric model. 

We also could not use the PMF technique with all of our data, because during sev- 
eral observing sessions the radio emissions from IM Peg were too weak (< 1 mJy) for the 
star to be detected in any single scan. In several other sessions, the star was detected only 
intermittently in single scans and in those cases only on the most sensitive baselines. Fur- 
thermore, our principal reference source, 3C 454.3, which we selected because of its close 
proximity on the sky to IM Peg (0.7° ang ular separ ation) and very high radio brightness, 



has a complex and evolving structure (see iPaper it ) that made necessary the inclusion of 



structure corrections in our model. The PMF technique cannot provide such corrections, 
so if we had relied upon the PMF technique for our astrometry, then we still would have 
needed to image 3C 454.3 separately and generate structure models for it. Likewise, IM Peg 



had h ighl y time- variable and sometimes complex structure (see Paper VII, iBietenholz et al. 



20 111 . and lLebach et al.lll999l ). so we would have had to image this source separately, too. 



Our solution to overcome the drawbacks associated with each of the PRM and PMF 
techniques was to merge the two approaches in a way that gave us the advantages and shed 
many of the disadvantages of both. 



2 In our implementation of the PRM technique, we used the data from only a single reference source, 
3C 454.3, for phase calibration. Thus we had no way to identify phase adjustments that had elevation- angle 
dependencies from those that did not. However, the PRM technique is not intrinsically limited to the use of 
one refere nce source for phase calib ration. In fact, software within the NRAO's AIPS package (AIPS task 
ATMCA: iFomalont fc Koganl 120051 1 is specifically intended to make use of multiple reference sources with 
the PRM technique. This software became available only at the end of our VLBI campaign, and due to time 
and budgetary constraints, we were unable to use it. The ability of this software to improve astr ometric 
accuracy, relative to when a single reference source is used, has been demonstrated ( Fomalont 20051) . 



5.2. Basic Description 



The basic approach behind the MA technique is first to follow the approach of the 
PMF technique to obtain model corrections. One bases these corrections on data from only 
those radio sources — usually the reference sources — whose emissions are sufficiently bright 
that they can be reliably detected in a single scan on many baselines and hence are well 
suited for analysis with the PMF technique. For improved accuracy, structure corrections 
for these sources can be obtained separately and included in the analysis. The data from 
the sources — usually the target sources — that are not well suited for the PMF technique are 
left unweighted but used as "placeholders" so that model corrections corresponding to the 
observation times and positions of those sources are readily attainable via interpolation of 
the model results obtained from the data that are weighted. In effect, more accurate models 
for the target-source observables are obtained than with the conventional PRM technique, 
because one has full use of the PMF-technique tools and models that are (in general) superior 
to those available with the PRM technique. The corrected models from the PMF technique 
are then used with the PRM technique, in place of the models and corrections that would 
have otherwise been used in a conventional implementation of the PRM technique. Our MA 
technique thus combines the superior model-correction capabilities of the PMF technique 
with the superior sensitivity of the PRM technique. 



5.3. Implementation 

Here we provide details about the way that we implemented the MA technique to analyze 
the IM Peg data obtained for GP-B. 



5.3.1. Observations 

Our observation strategy, as well as a list of the VLBI antennas we used, is provided in 



Paper III . We repeat some relevant points about our observations here. 



We had 35 sessions of VLBI observations of IM Peg between January 1997 and July 
2005. We made all of these observations over a continuum near 8.4 GHz (A = 3.6 cm). 
Typically we used data from 12 to 14 VLBI antennas located around the world. Most of the 
antennas we used are within the United States. 



We sequenced our observations through either three or four sources in a repeating 5.5- 
7 minute cycle. The three sources we observed in all 35 sessions were: our target source, 
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IM Peg; our principal extragalactic reference source, 3C 454.3 (B2251+158); and a secondary 
extragalactic reference source, B2250+194. We added the fourth source, B2252+172, to our 
cycle for the final 12 sessions. We treated this source as a second target source and used it 
to place bounds on the sizes of apparent changes in the positions of our reference sources 
due to structural evolution (i.e., changes in the brightness distributions of the reference 
sources). We also used it to test our different astrometric analysis techniques: We expected 
B2252+172 (unlike IM Peg) to be a stationary (extragalactic) object on the sky, so any 
estimated motions of the source would be a measure of the experimental errors inherent 
to our astrometric technique. Figure [1] shows the relative positions on the sky of the four 
sources we observed. Note that IM Peg is between, and approximately collinear with, the 
reference sources 3C 454.3 and B2250+194. We intentionally selected our reference sources 
with such alignment to obtain more accurate models for IM Peg data from the interpolation 
of the models fit to our reference-source data: These models account approximately for 
the elevation-angle dependencies of our measured fringe phases. We selected B2252+172 
primarily on the basis of its close proximity on the sky to 3C 454.3 and secondarily for its 
approximate alignment, too, between 3C 454.3 and B2250+194. 

The cycle time, and the scan time on each source per cycle, were chosen with the intent 
that they would be: (i) long enough to get reliable fringe-phase measurements for both 
reference sources on all baselines within a single scan; (ii) long enough to perhaps get reliable 
fringe-phase measurements of the guide star on at least the most sensitive baselines within a 
single scan, thereby allowing use of the PMF technique; and (iii) short enough to determine 
(without ambiguity) the integral number of 2tc phase wraps between successive measurements 
of each of the sources. We used two reference sources rather than one primarily so that 
we could distinguish model errors that have elevation-angle dependence (e.g., atmospheric 
delays) from those that do not (e.g., station clo ck behavior) and thereby improve our overall 
astrometric accuracy (see, too, 



Fomalontl 120051 1 . 



We also analyzed data from four sessions of VLBI observations of IM Peg made by one 
of us and other colleagues jLestrade et al.lll999T) betw een December 1991 and July 1994 



in 



support of the Hipparcos mission ( ILestrade et al.lll995f ). However, these four sets of observa- 



tions had significant differences from the later 35 sessions of observations made specifically 
for GP-B. Perhaps most significantly, these earlier observations used only four antennas per 
session and a single reference source, 3C 454.3. Therefore we did not include the results from 
these four earlier sessions in the comparisons of the astrometric techniques that we present 
in this paper. 
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B2252+1 72 



IM Peg 
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22.95 22.9 22.85 

Right Ascension (h) 



Fig. 1. — Positions (J2000) on the sky of the four radio sources used for GP-B astrometry. 
The east-west and north-south directions on the plot are shown to the same scale. 
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5.3.2. First-stage data reduction: calibration and fringe fitting 

We obtained the amplitude and phase calibration of our data for the MA technique by 
following steps 1-5 outlined in §[3] for the PRM technique. We then used AIPS task FRING 
to estimate the fringe phases, fringe rates, and group delays from the visibility (i.e., cross- 
correlation) data output from the correlator. We obtained these estimates on a scan-by-scan 
basis. We also obtained these estimates one baseline at a time by specifying two antennas 
at a time with the ANTENNA parameter of FRING Jf| The visibility data from which these 
estimates for each baseline were derived were independent, and thus we followed the common 
(and greatly simplifying) practice of treating the estimates for each baseline as independent. 

To account for the non-pointlike brightness distribution (i.e., structure) of our reference 
sources, we made i mages of the s ources using the self-calibration (or "hybrid mapping" ) 



scheme described by IWalkerl (119991 ) . We generated one such image for each reference source 
for each session of observations. We obtained structure corrections for our estimated fringe 
phases, fringe rates, and group delays by running AIPS task FRING twice for each baseline: 
once with the self-calibrated CLEAN map used as input calibration, and once with the source 
modeled as a point source (which is the FRING default). We took the structure corrections 
for the observables to be the differences between the estimates with a point-source map 
and with the CLEAN map. An important consideration in the generation of the structure 
corrections is the choice within the CLEAN map of the reference position that is defined to 
have zero structure correction (i.e., the effective point-source position of the source). Initially, 
we used the brightness peak of the image of each reference source as the reference position, 
because the brightne ss peak was always c learly defined and easily identifiable. However, as we 



discuss in Paper III (IBartel et al.l (120111 ); see also Paper II), further studies of the collection 
of images of our primary reference source 3C 454.3 revealed that a different point within 
the source made a substantially better choice of presumed- stationary reference position. We 
thus adopted this point, identified as "CI" in Papers II and III, as the reference position for 
our 3C 454.3 structure corrections. A set of images th at show th e location of "CI" within 



the brightness distribution of 3C 454.3 can be found in lPaper HI 



The results from AIPS task FRING for fringe phase, fringe rate, and group delay are 
estimates of residuals to the model used by the correlator. We used a customized version 
of AIPS task CL2HF0 to add these residual estimates to the a priori model values from the 



3 AIPS task BLING can also be used to estimate the observables one baseline at a time. However, technical 
difficulties with BLING in the earliest days of our AIPS processing of the experimental data compelled us 
to use FRING instead, and we then opted to continue to use FRING for all experimental data. 

4 Our customized version of this task performs the same basic functions. It includes features we required 
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correlator and place the resultant "total" observables into AIPS tables. We then "exported" 
these total observables from AIPS with AIPS task HF2SV. The resultant directory of files of 
total observables was combined with files of EOPs, surface meteorological data, antenna and 
source coordinates, and additional calibration data (see below) to produce databases that 
we used for the next stage of data reduction. We created these databases with the support 
of the VLBI group of the Space Geodesy Program at NASA's Goddard Space Flight Center. 
The theoretical observables and partial derivatives with respect to mode l parameters were 
computed with that VLBI group 's CALC software JCaprette etld1ll99(iL_We used CALC 
version 9.13 and the Jet Propulsion Laboratory's DE200 ephemeris ( jStandishl Il982l . 119901 ) 
for all of our a priori models. 

Along with the observable data, we also exported from AIPS a collection of files of 
the model fringe phases and rates used by the correlator. In the final stage of our data 
reduction, discussed in § 15.3.41 below, we effectively replaced these model phases and rates 
with improved models that we derived from the observable data that we exported. 



5.3.3. Second-stage data reduction: phase connection and parametric model fitting 

The next step of our data reduction was to use the measured phases from the reference 
sources to improve the model estimates of the phases of the target source. In particular, we 
used the measured phases from our two main reference sources, 3C 454.3 and B2250+194, 
to separate a priori model errors that have an elevation-angle dependence (as would be 
contained in r atm (u,t) in Equation [TJ from those that do not (as would be contained in 
Ti nst (a;,t) in Equation [[J. We also used these measured phases to identify time intervals 
when the phases cannot be well tracked by existing models. We unweighted the target- 
source data within these intervals, which could be as short as a single scan of target-source 
data. 

"Phase connection" refers to the determination of the correct integral number of 2n 
phase wraps between successive measurements of fringe phase for a particular baseline and 
source (see Figure [2]). Phase connection is possible in part because of our very accurate a 
priori models and in part because of the very stable time reference (hydrogen-maser frequency 
standard) used at each VLBI antenna. Our basic approach to phase connection was to 
start with the baselines and time segments of 3C 454.3 and B2250+194 data for which 
our a priori models were sufficiently accurate that the number of 2n phase wraps between 
successive fringe-phase measurements was immediately evident (see Figure [3]). We then used 



(and added ourselves) that were not available in the earliest versions of CL2HF. 
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this subset of phase-connected data to estimate adjustments to the parameters of our a priori 
model. The resulting improvement in our model enabled us to phase connect more data, 
and so the process proceeded in an iterative manner until we could no longer reliably phase 
connect any additional data. For observations above 10° elevation, we reliably connected, in 
total, 83% of the 3C 454.3 and B2250+194 phases that we evaluated in this data-reduction 
stage. 

The terms in Equation [T] that can be associated with difficulties in the phase-connection 
process are those for which the model errors can result in relatively large (> 1 radian) 
residual phase fluctuations in a random pattern over successive reference-source scans, i.e., 
in just a few minutes of time. Our a priori values for parameters such as site coordinates 
and EOPs were sufficiently accurate§ that they could not cause such rapid residual phase 
fluctuations. In fact, nearly all such phase fluctuations could be attributed to changes either 
in instrumental delays, r inst (u;,£), which include possible "clock jumps" in the frequency 
standards at the sites, or, most commonly, in propagation delays through the atmosphere, 
T a tm(w,£)- The data that we could not reliably phase connect tended to involve sites with 
large wet atmospheric delays (the St. Croix VLB A site was generally the most problematic 
for phase connection) or be from observations just above the 10° elevation threshold. 

We used surface meteorological data ( barometric pressure, temperature, and humidity 
or dew point) along with the equations from lSaastamoinenl (119721 ) to obtain a priori estimates 
°f Tzen_dry(0 and T 7iPn w P t(t) . For the map ping functions m dry (t,e) and m wet (t,e) we used the 
formulas provided by iNielll (119961 120001 ) . Our values for r mn (oj, t, s) came from the United 
States Air Force's Parameterized Ionos phere Model (P IM; iDaniell et al. 1995 ). which we 
adapted for use with VLBI observables JCampbelllll999l )FI 



In fitting our data to a model, we used a Kalman-filter estimator ( "SOLVK" ; iHerring et al. 



19901 ) that allowed us to model atmospheric delays at zenith and instrumental drifts as 
Gauss-Markov stochastic processes. Our VLBI data do not provide sufficient information to 



Our a priori values for antenna coordinates came primarily from the ITRF2000 solution (jAltamimi et al 



20021 ) provided by the International Earth Rotation Service (IERS). Our a priori values for the EOPs also 



all came from the IERS. 

6 One can obtain more accurate (we found) values for Ti on (oj,i,s) via the use of AIPS task TECOR and 
publicly available maps of to tal electron content (TEC) derived from Global Positioning System (GPS) data 



(jWalker fc Chatterieel _999). However, we found no useful TEC maps for our VLBI observations prior to 



September 1998. To avoid introducing possible systematic errors into our astrometric results from the use of 
different ionosphere models during different periods of observations, we opted to use PIM for all of our VLBI 
experiments. Ionosphere model (i.e., PIM) errors contribute insignificantly to th e total standard errors in 



the astrometric results for GP-B that we present in Paper V (jRatner et al.ll2011l) 
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Time (h UT on 2005 July 16) 

Fig. 2. — Representative phase residuals at the start of the phase-connection process. Shown 
are the differences between the measured fringe phases and the a priori model phases for the 
Pie Town (New Mexico) to North Liberty (Iowa) baseline of VLBI observations made on 2005 
July 16. Here the integral number of phase wraps for each measurement was selected such 
that the residual phase is always between 0° and 360° (i.e., between and 2n radians). We 
plot only the data for our two principal reference sources, 3C 454.3 and B2250+194. Error 
bars are not shown to improve plot clarity, but in general are <C 1 radian. The section of data 
with solid (for 3C 454.3) and dotted (for B2250+194) lines through the points corresponds 
to the period when changes in the model phase errors between successive measurements 
were <C 2n radians for both sources, so that the integral number of 2ir phase wraps between 
the measurements could be reliably determined. Only this "lined" section of the data is 
considered "phase connected" at this stage of the processing. 
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Fig. 3. — The same residual phases as shown in Figure |2] for the Pie Town (New Mexico) 
to North Liberty (Iowa) baseline after the integer number of 360° (or 2n radian) phase 
wraps were adjusted to remove the obvious 2tt radian phase jumps. The successive dashed 
horizontal lines are separated by 2n radians. The data are considered to be reliably phase 
connected when the changes in the residual phases between successive measurements are 
<C 2-7T radians. Only the data points with solid lines through them (one line for the 3C 454.3 
data and one line for the B2250+194 data) are assumed to be reliably phase connected at 
this stage. The phase-connected residuals still span a range > 2ix radians; this result is due 
to inaccuracies in the preceding model fit that a subsequent model fit would, in general, 
remove. The newly phase-connected data are used in a subsequent model fit to estimate 
corrections to the previous model and thereby facilitate the phase connection of additional 
data. This process can be iterated until no more 2ir phase wraps yield residual phase changes 
<C 27r radians between successive measurements. 
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simultaneously estimate adjustments to the separate terms in Equation [3j thus we adjusted 
only r zen _ wet (t), the term we assumed to have the least accurate a priori model. Figure 0] 
shows the adjustments to r zcn _ wet (t) that we obtained with our Kalman- filter estimator for 
the data corresponding to Figures |2] and [3j 

Through the use of data from two reference sources, 3C 454.3 and B2250+194, we were 
able to distinguish model adjustments to r atm (w,t), which has elevation-angle dependencies, 
from model adjustments to other terms, e.g., Ti nst (uj,t), which do not. Large phase jumps 
between successive reference-source scans were generally attributed to atmospheric fluctua- 
tions. Since the size of atmospheric fluctuations roughly scaled with the signal path length 
through the atmosphere, phase connection was generally easiest when the observed sources 
were at high elevation angles. Thus we generally started our phase- connect ion process with 
data obtained during the middle of an observing session, when the sources were at their 
highest elevation angles at antennas near the middle of the array. We then worked our way 
"outward" towards the early and late scans of the session. 

When we could not confidently determine the integral number of 2n phase wraps be- 
tween consecutive reference-source scans, we unweighted the target-source data between 
those scans. We also inserted a "break" marker at the point of the undetermined phase 
jump; the Kalman-filter estimator takes account of such breaks in its model fit. 

Figure |5]illustrates target-source residuals to the phase models derived from the reference- 
source data. The integral numbers of 2n phase wraps between successive scans are reliably 
determined for all weighted data. 



5.3.4- Final-stage data reduction: use of improved model in phase-referenced maps 



Once we had phase-delay and fringe-rate models based on parameter fits to the reference- 
source data, we imported those models into AIPS to use as the basis for our phase-referenced 
maps of the target source, either IM Peg or B2252+172. First we created files of fringe-phase 
and fringe-rate differences between the models derived from the reference-source data and 
the models (from the correlator) originally used in AIPS. Then we used AIPS task TBIN to 
create an AIPS "solution" (SN) table with these differences. We could then treat the change 
from the original models used in AIPS to the improved models derived from our reference- 
source data as a standard AIPS calibration step. We used AIPS task CLCAL to incorporate 
these model adjustments. We the n made phase-referenced maps with AIPS task IMAGR 



following standard procedures (e.g.. lBeasley fc Conwaylll995l ). We obtained positions for our 



target sources from these IMAGR maps as described in Paper VI (IRansom et al.l 12011m 
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Fig. 4. — Adjustments to model "wet" atmospheric delays at zenith (r zen _ wet ) estimated from 
the data with our Kalman filter for two VLBI sites, Pie Town (NM) and North Liberty (IA), 
for our 2005 July 16 observations. For scale, a 2n radian wrap in fringe phase corresponds to 
119 ps of delay along the line of sight from each antenna to each observed source. Although 
the adjustments are modeled as wet delays at zenith, they include contributions from model 
errors in the ionosphere, too, since those errors also have elevation-angle dependence. 
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Fig. 5. — Sample phase residuals for source B2252+172. The phase model was derived from 
our reference- source (3C 454.3 and B2250+194) data with our Kalman-filter estimator. The 
residuals shown are for the North Liberty (Iowa)-Pie Town (New Mexico) baseline and are 
from our 2005 July 16 observations. The short gaps occasionally seen between adjacent data 
points are due (in general) either to insufficient signal strength from B2252+172 to detect 
the source in a single scan of data, or to difficulties in determining the integral number of 
2n phase wraps between adjacent reference- source data. The number of 2n phase wraps 
chosen for each point was the number that minimized the magnitude of the residual. These 
magnitudes, and the differences in phase residuals between adjacent scans, are all <C 2n 
radians (<C 360°), which demonstrates that the data are properly phase connected. We used 
a relatively conservative estimate for the error of each phase measurement, so the scatter of 
the points relative to the size of the error bars appears somewhat small. 
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6. Comparison of Results 

We used two figures of merit to compare the results of our MA technique with the more 
conventional PRM and PMF techniques: the quality, as described below, of the B2252+172 
and IM Peg images produced (MA and PRM techniques only); and the level of consistency 
of the astrometric results for each of these two sources. 



6.1. Comparison of images 

Figure [6] shows a representative map of B2252+172 produced via the conventional PRM 
technique alongside the corresponding map generated via our MA technique. Table [1] lists 
the brightness-peak flux density, the dynamic range (i.e., the ratio between the flux density 
of the brightness peak and the flux density of the image noise floor), and the ratio of the 
peak-brightness amplitude (^4 p k+) to the amplitude of the most negative peak (A p k_) that we 
obtained with the PRM and MA techniques for each of our 12 sessions of observations of this 
source. If the extent of the brightness distribution of a source is comparable to, or larger than, 
the size of the restoring beam of its image, then the measured peak brightness and (hence) 
the dynamic range are functions of the size of the restoring beam. The size (and orientation) 
of the restoring beam can in turn depend upon which data were weighted (among other 
factors) when the image was produced. Since we employed different data-selection criteria 
with the PRM and MA techniques and hence did not necessarily use identical data with 
the two techniques, we consider the A pk+ /A pk _ ratio to be the most meaningful measure of 
image quality among the three Table [T] metrics. 

For both the PRM and MA techniques, we show these metrics for images produced after 
we applied either a 10° or a 30° elevation-angle cutoff to the data (i.e., after we unweighted all 
data from each site whenever the observed source was below the specified elevation angle). 
Since the major source of phase-calibration error with the PRM technique is often poor 
modeling of the atmosphere, and since (as mentioned earlier) atmospheric model errors tend 
to get larger at low elevation angles, the removal of low-elevation-angle data can potentially 
improve the quality of both the image and the astrometry obtained with the PRM technique. 
(Note that, for the PRM technique, the images when a 30° elevation-angle cutoff is used 
nearly always have higher dynamic ranges and A pk+ /Ap k ^ ratios than the corresponding 
images when a 10° elevation-angle cutoff is used.) For the MA technique, only those time 
spans with reliably phase-connected data are weighted, so we might expect a degradation of 
image quality (due to the reduction of weighted data) when a higher elevation-angle cutoff 
is imposed. The results found for our images are mostly consistent with this expectation, 
although Table [1] shows that for B2252+172, the higher elevation- angle cutoff did improve 
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Table 1: Characteristics of B2252+172 images obtained from the PRM and MA techniques. 



PRM Technique MA Technique 

With 10° With 30° With 10° With 30° 

elev. -angle cutoff elev. -angle cutoff elev. -angle cutoff elev. -angle cutoff 



Session 




Dyn. 


A pk + ° 




Dyn. 


A pk+ 




Dyn. 


A pk+ 




Dyn. 


A pk + 


Start Date 




Range b 


A pk- 


A pk+ 


Range 


A pk- 




Range 


A pk- 




Range 


A pk- 


2002-11-20 


12.2 


64 


9.7 


14.2 


68 


13. 


15.2 


165 


37. 


15.5 


170 


31. 


2003-01-26 


11.3 


53 


5.6 


14.1 


70 


8.3 


17.0 


140 


29. 


17.9 


149 


38. 


2003-05-18 


12.5 


62 


7.8 


16.6 


60 


8.7 


19.1 


187 


47. 


19.8 


173 


44. 


2003-09-08 


11.7 


64 


13. 


13.7 


72 


16. 


15.4 


141 


34. 


16.0 


132 


34. 


2003-12-05 


11.0 


51 


7.2 


14.3 


61 


9.9 


15.4 


106 


16. 


16.0 


97 


20. 


2004-03-06 


8.9 


56 


8.8 


11.0 


76 


12. 


12.1 


124 


30. 


12.4 


110 


25. 


2004-05-18 


6.5 


37 


7.3 


8.8 


48 


8.8 


10.0 


77 


18. 


10.5 


68 


18. 


2004-06-26 


7.0 


54 


11. 


9.0 


70 


14. 


9.8 


109 


25. 


10.1 


103 


22. 


2004-12-11 


9.0 


58 


8.9 


12.0 


94 


14. 


12.2 


128 


34. 


12.8 


116 


30. 


2005-01-15 


9.7 


87 


12. 


11.0 


103 


12. 


11.3 


175 


39. 


11.6 


149 


38. 


2005-05-28 


11.7 


48 


10.0 


15.5 


69 


15. 


16.3 


125 


36. 


16.5 


114 


33. 


2005-07-16 


12.5 


64 


11. 


15.4 


80 


14. 


16.0 


127 


27. 


16.5 


131 


28. 


MEAN d 


10.3 


58.1 


9.3 


13.0 


72.6 


12.0 


14.2 


133.6 


31.0 


14.6 


125.9 


30.0 



"Brightness-peak amplitude, in mJy/beam. 
b Dynamic range (see text). 

c Ratio of the brightness-peak amplitude (^4 p k+) to the amplitude of the most negative peak (A pk _). 
^Unweighted mean over all sessions. 
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Fig. 6. — Representative images of B2252+172 created in AIPS with the PRM technique 
(left) and the MA technique (right). The data for the images are from our observations on 
2004 March 6-7. Note that the two images have different coordinate origins. For each image, 
the size and orientation of the restoring beam is shown in the box in the lower left corner. 
The image created with the MA technique is visibly cleaner even with a lower minimum 
brightness contour. On average over all of our sessions of observations, the MA images had 
1.9 times the dynamic range of the PRM images. 
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the A p k + /A p k^ ratio with the MA technique in four of the twelve sessions. In any event, even 
if, for the PRM technique, we select for each session the elevation- angle cutoff in Table [1] 
with the highest dynamic range or A p k+ / A p ^_ ratio for each session, the corresponding MA 
images for a 10° elevation-angle cutoff have, on average, 1.9 times the dynamic range and 
2.7 times the A p k+/Apk- ratio of the PRM images. The MA images, when compared with 
the PRM images, also show significantly less variability in peak flux density for the two 
different elevation- angle cutoffs that are used. This superior consistency, too, suggests that 
the MA technique provides more accurate model calibration than the PRM technique. 

Figure [7] shows representative images of IM Peg produced with both the PRM and MA 
techniques. Averaged across all of our observing sessions, the peak flux density we detected 
for IM Peg was 16% higher with the MA technique than with the PRM technique. However, 
the significance of this result is unclear, because we used different restoring beams for each 
technique. Unfortunately, the data files used to make the images of IM Peg with the PRM 
technique were no longer available at the time of preparation of this publication, so we could 
not use restoring beams with the same size and orientation for both techniques. We also did 
not have noise-floor or "minimum peak" information available for these PRM images. In 
addition, we know that IM Peg can vary in total brightness and brightness distribution (i.e ., 
radio structure) over hour time scales (Paper VII. iBietenholz et al.ll2011uLebach et al.lll999l ). 
which further complicates comparisons between the techniques, especially since a constant 
flux density over the full duration of a session of observations is inherently assumed in the 
imaging process for both techniques. Thus, for comparison of image quality from the PRM 
and MA techniques, we considered the B2252+172 results a better measure than the IM Peg 
results. 



6.2. Comparison of astrometric results 

Figure |8] compares the position estimates for B2252+172 (relative to 3C 454.3) from the 
PRM, PMF, and MA techniques. (We used a 30°, rather than 10°, elevation-angle cutoff 
of data with the PRM technique, because the higher cutoff yielded better results for this 
technique. We used our nominal 10° elevation- angle cutoff of data with the PMF and MA 
techniques.) Table [2] shows the root-mean-square (rms) scatter about the mean position for 
each technique. When, as with B2252+172, the target source is commonly detected on many 
baselines in a single scan and structure corrections with the PMF technique are available, the 
MA technique shows no astrometric advantage over the PMF technique. However, our MA 
technique yields position estimates with only about 60% of the rms scatter of the estimates 
from the PRM technique, i.e., removes the equivalent of a noise source of over 0.06 mas per 
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Fig. 7. — Representative images of IM Peg created in AIPS with the PRM technique (left) 
and the MA technique (right), from observations on 2004 March 6-7. For each image, the 
size and orientation of the restoring beam are shown in the box in the lower left corner. Note 
that the left and right images have different coordinate origins. 
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coordinate estimate associated with the PRM technique. 

The model for IM Peg's position over time must include proper motion, parallax, and 
orbital moti ons associated with the binary system. We used the preferred parametrization 



of Paper V (IRatner et al.ll201ll ). one in which no proper acceleration is estimated, to obtain 



the results presented here. 

For seven of the 35 observing sessions of IM Peg, including the most astrometrically 
important final three, the radio emissions from the star were insufficiently bright to obtain 
reliable detections of IM Peg in single scans. Thus we could not obtain position estimates 
of IM Peg with the PMF technique for any of these sessions. We also found from images 
obtained with the PRM and MA techniques that IM Pe g had relatively complex structure 



during several of our observing sessions (see iPaper VIII ). and the PMF technique by itself 



cannot well identify or characterize such structure. Thus we deemed the PMF technique 
unsuitable for astrometric studies of IM Peg with our data. 

In Figure [9] we show the post-fit residuals of the position estimates for IM Peg from the 
PRM and MA techniques, and in Table [3] the rms scatter of these same pos-tfit residuals. 
On average the post-fit residuals from our MA technique had about 10% less scatter in right 
ascension and 7% less scatter in declination than the corresponding residuals from the PRM 
technique. The scatter from both techniques is likely dominated by the unmodel ed intrinsi c 



motions of the stellar radio emission relative to the star's primary component ( IPaper VII ). 
Hence, since we have no model that can accurately describe the motion of the radio emission 
of IM Peg relative to its primary, we can expect at most only a minor improvement in the 
rms scatter of position residuals from the use of any superior astrometric technique. Indeed 
the MA technique yielded just such a minor improvement. 



7. Conclusions 

In comparing the astrometric results from phase-referenced mapping (PRM), parametric 
model fitting (PMF), and the merged analysis (MA) technique introduced in this paper, we 
found: 

1. The PRM technique can be used with relatively weak radio sources. It is the least 
labor intensive among the techniques tested, but it also provides the least astrometric 
accuracy. (With this technique, unlike with the PMF or MA techniques, we were unable 
to use multiple reference sources to better account for the elevation-angle dependence 
of model errors; such use could improve the astrometric accuracy attainable with this 
technique.) 
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Fig. 8. — Position estimates for B2252+172 relative to 3C 454.3 from the PRM (diagonal 
crosses; connected by dotted line segments), PMF (diamonds; connected by dashed line 
segments), and MA (filled circles; connected by solid line segments) techniques (see text). 
The zero position for each coordinate is the unweighted mean value of the twelve estimates 
obtained with the MA technique. For better plot clarity, the error bars for the plotted 
points are not shown. Since the scatter in the position estimates is due mainly to model 
errors (such as for the ionosphere) rather than to SNR limits of the data, the rms scatter of 
the position estimates relative to the mean for each technique, as given in Table [21 is a good 
approximation for the errors of the corresponding points shown in this plot. Such errors 
ignore possible systematic effects that could have shifted the mean positions for each of the 
three techniques. 
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Table 2: Astrometric results for B2252+172 from the PRM, PMF, and MA techniques. 





Mean Offset a 


RMS Scatter 15 


Standard Error 




(mas) 


(mas) 


of mean c (mas) 


Technique 


a 5 


a S 


a 


S 


PRM 


-0.039 0.021 


0.082 0.092 


0.025 


0.028 


PMF 


0.012 0.019 


0.048 0.053 


0.014 


0.016 


MA 


0.000 0.000 


0.050 0.057 


0.014 


0.017 



a Dcfined to be zero for the MA technique in both right ascension (a) and declination (S). The mean offsets 
for the PRM and PMF techniques are relative to the mean coordinate position from the MA technique. 
^Relative to the mean coordinate position for each technique. 
Calculated from the rms scatter. 

Note. — For each technique we assumed, for each of a and 5, the same standard error for every observing 
session. The rms scatter for the PRM technique is the largest; the levels of rms scatter for the PMF and MA 
techniques are not significantly different from each other, although the former appear to be slightly smaller. 



Table 3: RMS scatter of IM Peg position residuals obtained from the PRM and MA tech- 
niques. 





RMS Scatter 




(mas) 


Technique 


a 5 


PRM 


0.40 0.44 


MA 


0.35 0.41 



Note. — The position residuals are relative to an astrometric model fit in which the radio emissions are 
assumed to be tied to a single stellar component of the IM Peg binary system. We believe that the rms 
scatter is due predominantly to intrinsic motions of the radio emissions relative to this stellar component 
rather than to measurement accuracy. Nevertheless, the rms scatter is roughly 10% smaller with the MA 
technique than with the PRM technique. 
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Fig. 9. — Postfit residuals of the position estimates of IM Peg from the PRM (diagonal 
crosses; connected by dotted line segments) and MA (filled circles; connected by solid line 
segments) techniques. (In the model fit, we estimate the position, proper motion, parallax, 
and apparent binary-system orbit of the stellar radio emitting region.) We could not obtain 
position estimates of IM Peg with the PMF technique for seven of our 35 observing sessions 
due to insufficient radio brightness of the star; hence we deemed the standard PMF technique 
unsuitable for astrometric studies of IM Peg with our data. For better plot clarity, the error 
bars for the plotted points are not shown. The scatter in the position estimates from both 
the PRM and MA techniques is due mainly to intrinsic motions of the radio emission of 
IM Peg relative to the primary component of the star rather than to errors in the position 
estimates of the radio emission. 



-30 - 



2. The PMF technique is useful only when both reference and target sources are: (i) re- 
liably detectable in short (i.e., few- minute) scans of data; and (ii) pointlike, unless 
structure corrections are available from other means. Under these circumstances, the 
PMF technique yielded the highest astrometric accuracy among the techniques tested. 

3. Our MA technique can provide astrometric accuracies nearly equal to those obtained 
with the PMF technique. 

4. Very importantly, our MA technique can be used with sources that are too weak for use 
of the PMF technique, as was IM Peg during seven of our 35 VLBI sessions. In general, 
any source that is sufficiently bright for use of the PRM technique is sufficiently bright 
for use of our MA technique. The MA technique can also be used with sources that 
have complex brightness distributions and was shown to yield images with significantly 
higher dynamic ranges than comparable images from the PRM technique. 

Our use of the MA technique enabled us to obtain more accurate astrometric measurements 
of IM Peg, the guide star for the GP-B mission, than we could have otherwise obtained with 
the conventional PRM or PMF techniques. 

ACKNOWLEDGMENTS. We thank N. Nunes for developing the graphical software 
package we used to phase connect our VLBI data. We thank the VLBI group of NASA's 
Goddard Space Flight Center, L. Petrov in particular, for providing technical support and 
information. We are also grateful to the many people involved in our campaign of VLBI ob- 
servations who went above and beyond the call of duty, including S. Dains, C. Garcia Miro, 
E. Moll, and L. Cameron. We thank C. Jacobs and O. Sovers for information about the coor- 
dinates of the three 70 m NASA DSN antennas used in our observations, and R. C. Walker for 
information about the VLA coordinates. We thank J. L. Davis and S. S. Shapiro for their in- 
formation and insights on the use of our Kalman-filter estimator (SOLVK). Finally, we thank 
R. C. Walker, K. Desai, and E. Greisen for their support in our use of AIPS. The National 
Radio Astronomy Observatory (NRAO) is a facility of the National Science Foundation op- 
erated under cooperative agreement by Associated Universities, Inc. This research made use 
of NASA's Astrophysics Data System, which was conceived, developed, and continues to be 
operated by the Smithsonian Astrophysical Observatory at the Harvard- Smithsonian Center 
for Astrophysics. Our work was supported by NASA prime award NAS8-39225, Stanford 
University sub-award PR 6750, the Smithsonian Institution, and York University. 



31 



REFERENCES 

Altamimi, Z., Sillard, P., & Boucher, C. 2002, J. Geophys. Res. (Solid Earth), 107, 2 

Bartel, N., Bietenholz, M. F., Lebach, D. E., Lederman, J. I., Petrov, L., Ransom, R. R., 
Ratner, M. I., & Shapiro, I. I. 2012, this issue (Paper III) 

Bartel, N., Ratner, M. I., Shapiro, I. I., Cappallo, R. J., Rogers, A. E. E., & Whitney, A. R. 
1985, AJ, 90, 318 

Beasley, A. J., & Conway, J. E. 1995, in ASP Conf. Ser. 82: Very Long Baseline Interferom- 
etry and the VLBA, ed. J. A. Zensus, P. J. Diamond, & P. J. Napier, 328 

Bietenholz, M. F., Bartel, N., Lebach, D. E., Ransom, R. R., Ratner, M. I., & Shapiro, I. I. 
2012, this issue (Paper VII) 

Campbell, R. M. 1999, New A Rev., 43, 617 

Caprette, D. S., Ma, C, & Ryan, J. W. 1990, NASA Technical Memorandum 100765 (Crustal 
Dynamics Project Data Analysis, 1990), National Aeronautics and Space Adminis- 
tration, Goddard Space Flight Center, Greenbelt, MD, Tech. rep. 

Clark, T. A., et al. 1985, IEEE Transactions on Geoscience and Remote Sensing, 23, 438 

Counselman, III, C. C, Kent, S. M., Knight, C. A., Shapiro, I. I., Clark, T. A., Hinteregger, 
H. F., Rogers, A. E. E., & Whitney, A. R. 1974, Phys. Rev. Lett., 33, 1621 

Daniell, R. E., Brown, L. D., Anderson, D. N., Fox, M. W., Doherty, P. FL, Decker, D. T., 
Sojka, J. J., & Schunk, R. W. 1995, Radio Science, 30, 1499 

Davis, J. L., Herring, T. A., Shapiro, I. I., Rogers, A. E. E., & Elgered, G. 1985, Radio 
Science, 20, 1593 

Diamond, P. J. 1995, in ASP Conf. Ser. 82: Very Long Baseline Interferometry and the 
VLBA, ed. J. A. Zensus, P. J. Diamond, & P. J. Napier, 227 

Dodson, R., Legge, D., Reynolds, J. E., & McCulloch, P. M. 2003, ApJ, 596, 1137 

Fey, A. L., Gordon, D., & Jacobs, C. S., (eds.) 2009, IERS Technical Note No. 35, 
|http://www.iers.org/IERS/EN/Publications/TechnicalNotes/tn35.html ( Frankfurt : 
Verlag des Bundesamts fur Kartographie und Geodasie) 

Fomalont, E. B. 2005, in ASP Conf. Ser. 340: Future Directions in High Resolution Astron- 
omy, ed. J. Romney & M. Reid, 460 



-32 - 



Fomalont, E. B., & Kogan, L. 2005, MPS Memo 111 (ATMCA: 
Phase Referencing Using More Than One Calibrator), Na- 
tional Radio Astronomy Observatory, Socorro, NM, Tech. rep., 
ftp77ftp.aoc.nrao.edu/pub/software/aips/TEXT/PUBL/AIPSMEMlll.PS| 

Gorenstein, M. V., et al. 1983, Science, 219, 54 

Herring, T. A., Davis, J. L., & Shapiro, I. I. 1990, J. Geophys. Res., 95, 12561 
Herrnstein, J. R., et al. 1999, Nature, 400, 539 

Ishihara, Y., Nakai, N., Iyomoto, N., Makishima, K., Diamond, P., & Hall, P. 2001, PASJ, 
53, 215 

Kondratko, P. T., Greenhill, L. J., & Moran, J. M. 2004, in Proc. 222nd IAU Symp., ed. 
T. Storchi-Bergmann, L. C. Ho, & H. R. Schmitt, 325-326 

Lebach, D. E., Ratner, M. I., Shapiro, I. I., Ransom, R. R., Bietenholz, M. F., Bartel, N., & 
Lestrade, J.-F. 1999, ApJ, 517, L43 

Lestrade, J.-F., Preston, R. A., Jones, D. L., Phillips, R. B., Rogers, A. E. E., Titus, M. A., 
Rioja, M. J., & Gabuzda, D. C. 1999, A&A, 344, 1014 

Lestrade, J.-F., Rogers, A. E. E., Whitney, A. R., Niell, A. E., Phillips, R. B., & Preston, 
R. A. 1990, AJ, 99, 1663 

Lestrade, J.-F., et al. 1995, A&A, 304, 182 

Ma, C, et al. 1998, AJ, 116, 516 

Mathews, P. M., & Shapiro, I. I. 1995, Physics of the Earth and Planetary Interiors, 92, 99 

Niell, A. E. 1996, J. Geophys. Res., 101, 3227 

— . 2000, Earth, Planets, and Space, 52, 699 

Pradel, N., Chariot, P., & Lestrade, J.-F. 2006, A&A, 452, 1099 

Ransom, R. R., Bartel, N., Bietenholz, M. F., Lebach, D. E., Lederman, J. I., Luca, P., 
Ratner, M. I., & Shapiro, I. I. 2012a, this issue (Paper II) 

Ransom, R. R., Bartel, N., Bietenholz, M. F., Lebach, D. E., Lestrade, J.-F., Ratner, M. I., 
& Shapiro, I. I. 2012b, this issue (Paper VI) 



-33 - 



Ratner, M. I., Bartel, N., Bietenholz, M. F., Lebach, D. E., Lestrade, J.-F., Ransom, R. R., 
& Shapiro, I. I. 2012, this issue (Paper V) 

Rodin, A., & Sekido, M. 2002, in Proc. 6th EVN Symp., ed. E. Ros, R. W. Porcas, A. P. 
Lobanov, & J. A. Zensus, 247 

Rogers, A. E. E., et al. 1974, ApJ, 193, 293 

Ryan, J. W., Ma, C, & Caprette, D. S. 1993, NASA Technical Memorandum 104572 (Data 
analysis - 1992, Final Report of the Crustal Dynamics Project VLBI Geodetic Results, 
1979 - 1991), National Aeronautics and Space Administration, Goddard Space Flight 
Center, Greenbelt, MD, Tech. rep. 

Saastamoinen, J. 1972, in The Use of Artificial Satellites for Geodesy, ed. S. W. Henriksen, 
A. Mancini, & B. H. Chovitz, 247 

Seidelmann, P. K. 1982, Celestial Mechanics, 27, 79 

— . 1992, Explanatory Supplement to the Astronomical Almanac (University Science Books) 

Shapiro, I. I., Bartel, N., Bietenholz, M. F., Lebach, D. E., Lestrade, J.-F., Ransom, R. R., 
& Ratner, M. I. 2012, this issue, (Paper I) 

Shapiro, I. I., et al. 1979, AJ, 84, 1459 

Shapiro, S. S., Davis, J. L., Lebach, D. E., & Gregory, J. S. 2004, Phys. Rev. Lett., 92, 
121101 

Standish, Jr., E. M. 1982, A&A, 114, 297 
-. 1990, A&A, 233, 252 

Walker, R. C. 1999, in ASP Conf. Ser. 180: Synthesis Imaging in Radio Astronomy II, ed. 
G. B. Taylor, C. L. Carilli, & R. A. Perley, 433 

Walker, R. C, & Chatterjee, S. 1999, VLBA Scientific Memo 23 (Ionospheric Corrections 
Using GPS Based Models), National Radio Astronomy Observatory, Socorro, NM, 
Tech. rep., |http:/ /www.vlba.nrao.edu/memos/sci/] 



This preprint was prepared with the A AS IATgX macros v5.2. 



